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Abstract 

We propose an ab initio method to evaluate the core-valence-valence (CVV) Auger spectrum of systems 
with filled valence bands. The method is based on the Cini-Sawatzky theory, and aims at estimating the 
parameters by first-principles calculations in the framework of density-functional theory (DFT). Photoe- 
mission energies and the interaction energy for the two holes in the final state are evaluated by performing 
DFT simulations for the system with varied population of electronic levels. Transition matrix elements are 
taken from atomic results. The approach takes into account the non-sphericity of the density of states of 
the emitting atom, spin-orbit interaction in core and valence, and non quadratic terms in the total energy 
expansion with respect to fractional occupation numbers. It is tested on two benchmark systems, Zn and Cu 
metals, leading in both cases to L23M45M45 Auger peaks within 2 eV from the experimental ones. Detailed 
analysis is presented on the relative weight of the various contributions considered in our method, providing 
the basis for future development. Especially problematic is the evaluation of the hole-hole interaction for 
systems with broad valence bands: our method underestimates its value in Cu, while we obtain excellent 
results for this quantity in Zn. 

PACS numbers: 79.20.Fv, 71.15.Mb, 82.80.-d 
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I. INTRODUCTION 



The Auger spectroscopy involves the creation of two localized holes at or close to the same 
atom, hence giving access to local electronic properties. Direct information on the local density 
of valence states is brought by core-core- valence transitions; in the case of core-valence- valence 
(CVV) ones, which will be investigated here, one can in addition access the screened Coulomb 
repulsion amongst the two valence holes in the final state which is relevant to a wide class of 
phenomena, and study its effects. 

From the theoretical point of view, a large amount of work has been devoted to the calculation 
of the Auger spectra of solids during the last three decades. - A general formulation of the dynam- 
ical Auger decay, where the creation of the initial core hole and the Auger decay are considered as 
coherent processes, was given by Gunnarsson and Schonhammer^ but is of hard practical imple- 
mentation. In solids with (almost) closed valence bands, where no dynamical core-hole screening 
can occur before the Auger decay, one can employ a simpler two-step approximation and consider 
the above events as independent. Under this assumption, Cini^ and Sawatzky'^ (CS) proposed a 
simple model providing the Green's function describing the two valence holes left after the Auger 
decay. Good agreement with experiments was achieved using fitting parameters for the screened 
Coulomb interaction, giving a quantitative understanding of the Auger spectra of transition metals 
located at the beginning and the end of the row, such as Ti,^ Ag,^ and Au.^ These results con- 
firmed the usefulness of including explicitly on-site Hubbard terms to one-body Hamiltonian, and 
prompted an extension to nearest- neighbour interactions .- 

These studies determined the relevant physical parameters by reproducing experimental find- 
ings within a semi-empirical approach. Of particular interest is the parameter governing the in- 
teraction amongst the two holes in the final state, which has an analogue in the popular LDA-i-U 
description for correlated systems.-"^ Even if methods for its ab initio evaluation have been pro- 
posed, also this quantity is often determined by phenomenological arguments. 

The possibility of evaluating CVV spectra from first-principles, rather than from a model with 
parameters fitted to experiments, would then be very desirable as it would allow predicting differ- 
ent situations (e.g., investigate the effect of a given chemical environment on the Auger current) 
and a deeper interpretation of experimental findings. 

This paper addresses such possibility, by proposing a method to compute the parameters en- 
tering the CS model by ab initio simulations. In this step towards a first-principle description of 
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CVV Auger spectra in systems where the interaction of the final-state holes cannot be neglected, 
we aim at highlighting the most important contributions to the spectrum, which one should fo- 
cus to in forthcoming improvements. The method is based on Density Functional Theory (DFT) 
simulations in the Kohn-Sham (KS) framework, with constrained occupations. We make use of 
comparison with reference atomic calculations to extrapolate the electronic properties of the sam- 
ple when they are more difficult to evaluate directly. Results are presented for the L23M45M45 
Auger line of Cu and Zn metals, which have been chosen as benchmark systems with closed 3d 
bands: the former being more challenging for the proposed procedure, and the latter bearing more 
resemblance with the atomic case. 

The paper is organized as follows. In section |ll] we describe our method to evaluate the Auger 
spectra by first-principles calculations. Section [ni] presents our theoretical results for Cu and Zn 
metals, comparing them with experimental results in the literature. In section |W] we analyze 
the weight of various contributions and discuss improvements. Finally, section |V] is devoted to 
conclusions. 

II. THEORETICAL METHODS 

A. Model Hamiltonian and the Cini-Savatzky solution 

We describe the electron system in the hole representation, by an Hubbard-like^ model Hamil- 
tonian 

^ = GccJCf -|-^£vcjcv + - ^ ^cpi9293(P4''Cpi'^92*'^4''93' 
V ^ <Pl92<P394 

where c and v label the core state involved in the transition and the valence states of the system, 
respectively, including the spin quantum number. In bulk materials, v is a continuous index. The 
last term is the hole-hole interaction Hamiltonian, parametrized by the screened repulsion U, and 
is for simplicity restricted here to a finite set of wavefunctions, cp, centered at the emitting atom 
(hence neglecting interatomic interactions). In closed shell systems, £c and yield the core and 
valence photoemission energies, since the two-body term has no contribution on the one-hole final 
state or on the zero-hole initial one. 

A two-step model is adopted to represent the Auger process, assuming that the initial ionization 
and the following Auger decay of the core hole can be treated as two independent events. In other 
terms, we assume that the Auger transition we are interested in follows a fully relaxed ionization 
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of a core shell. If the ground state energy of the neutral A'^-electron system is chosen as a reference, 
the energy of the initial state is simply given by £c. The total spectrum for electrons emitted with 
kinetic energy co is proportional to 

5(a))=£A^Dxy(£..-(0)Ay, (2) 

XY 

where X and Y are the final-state quantum numbers. Ax is the Auger matrix element corresponding 
to the final state X, and Dxy represents the two-hole density of states. Notice that Eq. ^ coincides 
with the Fermi golden rule if the states X, Y are eigenstates of the Hamiltonian, so that Dxy 
is diagonal. The presence of the transition matrix elements effectively reduces the set of states 
contributing to Eq. ^ to those with a significant weight close to the emitting atom. This motivates 
the approximation to restrict X and Y to two-hole states based on wavefunctions centered at the 
emitting atom, such as the set {cp} previously introduced. Therefore, the CVV spectrum is a 
measure of the two-hole local density of states (2hLD0S), with modifications due to the matrix 
elements. The 2hLD0S could in principle be determined as the imaginary part of the two-hole 
Green's function, Gxy, solution of Eq. ([T]). However, because of the presence of the hole-hole 
interaction term, evaluating Gxy is in general a formidable task. 

For systems with filled valence bands, the two holes are created in a no-hole vacuum and one 
is left with a two-body problem. A solution in this special case has been proposed by Cini^^ 
and Sawatzky,- and is briefly reviewed here (see Ref. [ij for an extended review). The two-holes 
interacting Green's function, G, is found as the solution to a Dyson equation with kernel U, which 
reads: 

G((o) = g(0)((o)(1-[/g(0)((d))"\ (3) 

Here, G^^) is the non-interacting Green's function which can be computed from the non-interacting 
2hLDOS,D(0), via Hilbert transform. Such 2hLD0S results from the self-convolution of the one- 
hole local density of states (IhLDOS), D^^^=d*d. 

The quantum numbers LSJMj (intermediate coupling scheme) are the most convenient choice 
to label the two-hole states, allowing for the straightforward inclusion of the spin-orbit interaction 
in the final state by adding to the Hamiltonian the usual diagonal term, proportional to [7(7+ 1) — 
L(L-M) - 5(5+ 1)]. Finally, the CVV lineshape is: 

G(o)(£,-a)) 



LSJMj 



l's'j'm' 



1-[/G(o)(£c-(o) 



ji r 



(4) 

LSJMj 

l's'j'm' I 



For comparison with experimental results, this is to be convoluted with a Voigt profile to account 
for core-hole lifetime and experimental resolution. 

It is customary to isolate two limiting regimes: (i) When U is small with respect to the valence 
band width W (broad, band-like spectra) the 2hLD0S is well represented by D^^\()i). However, 
in such a case it might be even qualitatively important to account for a dependence of the matrix 
elements on the Auger energy (O. As a consequence, accurate calculations of the lineshape require 
the simultaneous evaluation of the matrix elements and the DOS. (ii) For U larger than W , narrow 
atomic-like peaks dominate the spectrum, each peak from an LSJ component. Hence, to the first 
approximation the spectrum is described by a sum of 5-functions, weighted by matrix elements 
whose dependence on the Auger energy may be neglected. If we take the matrix U diagonal in the 
LSJ representation, and indicate by E^^^j the weighted average of Dj^^y((o), one obtains: 

5((0) ^ £(27+ l)\ALSj\'K^^c-Efl-ULSj)-^). (5) 

LSJ 

Atomic matrix elements can be taken as a first approximation, often satisfactory, and can be eval- 



uated as shown in Ref. 



I2I An approach which could bridge between these two limiting regimes. 



considering both finite values of U and the energy dependence of the matrix elements, is still 
missing to our knowledge. 

In the present work we adopt Eq. (HI) in order to simulate the spectrum. Accordingly, one has 
to determine the quantities A, U , D'^^\(£>), and e^. In this paper we make use of a i7 matrix which 
does not include the spin-orbit interaction, and is diagonal on the LS basis. We take atomic results 
in the literature for the matrix elements A, which are assumed independent of J too}^ The other 
quantities are computed by DFT simulations, as detailed in the following Section. 



B. Ab initio determination of the relevant parameters 

To evaluate ab initio the photoemission energies we use a method closely related to Slater's 
transition-state theory, while the parameter U is computed following a general procedure first 



proposed in Ref. 
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and then adopted by several authors. 
One extrapolates total energies for the system with N, N — I and N — 2 electrons by DFT 
calculations with constrained occupations for N — q electrons, with q small (typically, up to 0.05), 
so that ionized atoms in otherwise periodic systems can be treated in rather small supercells. We 
make the approximation that the total energy of the system with qt electrons removed from the 
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level / is given by a power expansion in qi up to third order: 

E{N-qi) =E{N) +A,q,+B,qj + Qql (6) 

In the following we shall assume that this can be extended to finite values of qi. The introduction 
of the cubic term Q^? allows for a ^-dependence of the screening properties of the system. The 
coefficients A/, 5,, and C/, where / labels core and valence states involved in the transition, are in 
this framework all what is needed to compute the Auger electron energy. They can be evaluated in 
two equivalent ways, whichever is most convenient: by taking the first, second, and third deriva- 
tives of the total energy E{N — qi) for qi — > 0; by using Janak's theoren>i^ and computing the KS 
eigenvalue of level i and its first and second derivatives: 

- eP (A^ -q,)= Ai -f 2Bm + 3Qqj. (7) 

In particular. A, is given by (minus) the KS eigenvalue in the neutral system. 

The binding energy of a photoemitted electron, ei=E]^^^ — E{N — 1,) — E{N), to be used in 
Eq. (dJ, is given by Eq. ^ as: 

£,=Ai + B,+Ci. (8) 

This is very close to the well-known Slater's transition- state approach, in which the XPS energy 
equals the (minus) eigenvalue at half filling. The latter amounts to A,- + Bi + ^Ci when approx- 
imating the total energy by a cubic expansion as in Eq. In other terms, it differs from the 
result of Eq. ([8]) only by |C,, with C, < 1 eV in the cases considered here (see below). It is worth 
noticing that the term 5, + Ci acts like a correction to the (minus) KS eigenvalue A,-, accounting 
for dynamical relaxation effects even though all terms are evaluated within KS-DFT. 

The evaluation of the A, B, and C coefficients for localized states poses no additional difficulty. 
Instead, care must be taken when determining those corresponding to the delocalized valence 
shells of bulk materials (Ay, By, and Cy), for which we propose the following method. As for 
Ay, this is a continuous function of the quantum number v and, by taking advantage of Janak's 
theorem, it is the KS band energy with reversed sign. To estimate By and Cy, we neglect their 
dependence on v and assume that a single value can be taken across the valence band, acting as a 
rigid shift of the band. Hence, the IhLDOS is obtained from the KS LDOS, (i^^((o), as: 

J((0) = J^^(-(0 + 5y+Cy). (9) 

For sake of the forthcoming discussion, one can also define a single value of Ay in the solid by 
taking the KS valence band average. 
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Is 2s 2p 3s 3p 


average 


Cu 


AS 
AC 


-4.77 -4.92 -4.90 -4.81 -4.76 
-0.88 -0.80 -0.82 -0.73 -0.72 


-4.85 ±0.07 
-0.81 ±0.07 


Zn 


AB 
AC 


-4.21 -4.26 -4.27 -4.19 -4.17 
-0.49 -0.33 -0.26 -0.32 -0.35 


-4.23 ±0.04 
-0.35 ±0.09 



TABLE I: Differences amongst the values of B and C in the bulk and the atom, AB = B — B" and AC = 
C — C, for core levels of Cu and Zn. The last column reports the average and standard deviation across the 
core levels. All values in eV. 

We expect the above approximation to be a good one as long as the valence band is sufficiently 
narrow and deep (since eventually the correction should approach zero at the Fermi level). Still 
under this simplification, the direct evaluation of By and Cy would ask for constraining the occu- 
pations for fairly delocalized states, which is a feasible but uneasy task. As an alternative route, 
we suggest a simpler approach based on the working hypothesis that the environment contribution 
to the screening of the positive charge qt in Eq. Q does not depend strongly on the shape of the 
charge distribution. Practically, we take the neutral isolated atom as a reference configuration in 
Eq. Q, and evaluate the coefficients 5" and for this system. The two quantities A5 = 5, — 
and AC = C, — Cf can be easily computed for core levels. Such bulk-atom corrections are re- 
ported in Table H] for Cu and Zn, which demonstrates that they are almost independent of the core 
level. This supports our working hypothesis, and enables us to extrapolate to the valence shell. 
Accordingly, and Cy are given by: 

B, = B^^+AB, (10) 
Cv = Q" + AC. (11) 

We remark that by choosing the neutral atom as the reference system some degree of arbitrari- 
ness is introduced. In principle, one could evaluate the atomic coefficients in a configuration which 
is closest to the one of the atom in the solid, depending on its chemical environment. However, 
such arbitrariness has limited effect on the final value of By (similar discussion applies for Cy), 
owing to cancellations between B" and B" in Eq. (flOl ). as will be demonstrated in the following. 

Regarding the interaction energy U for the two holes in a valence level, defined by [£'(A'^ — 
2) — ^(A^)] —2[E{N— 1) — ^(A^)], let us consider the case of spherically symmetric holes (non- 
spherical contributions, giving rise to multiple! splitting, will then be added). Such spherical 
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interaction, denoted by C/gph, can be determined via Eq. resulting in 

[/sph = 25v + 6Cv. (12) 

This amounts to the second derivative of the DFT energy as a function of the band occupa- 
tion, U{q) = d^E{N — q)/dq^, as originally suggested by Gunnarsson and coworkers, here 
evaluated for the — 1 -electron system rather than for the neutral one. Differently, the in- 
teraction energy commonly used in LDA-i-U calculations of the ground state is defined as 
E{N + 1) + E{N — 1) — 2E{N) and hence evaluated by the second derivative in q = 0, result- 
ing in 2By only. Notice here that the role of the cubic term in Eq. Q is to introduce a dependence 
of the interaction energy on the particle number, following the one of the screening properties of 
the system. Finally, non-spherical contributions, which give rise to multiplet splitting, are added to 
t/sph- It has been demonstrated^ for a number of materials, including Cu and Zn, that these terms 
are well reproduced by a sum of atomic Slater integrals, a2F^ + a^F^, where the coefficients 
and a4 depend on the multiplet configuration and 

noo poo yk 

F^= rldr rldr2^,[^^{n)<^%r2)]\ (13) 

Here (p'^(r) is the atomic radial wave function relevant to the process under investigation (e.g., the 
3d one for a CM^'-^M^s Auger transition), and r< (r>) is the smaller (larger) of ri and r2- Notice 
that the spherical Slater integral F^ is implicit into i/sphj which has the meaning of a screened 
Coulomb integral- 
Summarizing, one has: 

C/ = 25^ + 6C^f-F2A5 + 6AC + a2F2 + a4F4. (14) 

It is customary to write U = F — R, where F = F° -|- ajF^ -f a^F^, and R is the "relaxation 
energy". This can be further decomposed into an atomic and an extra-atomic contribution, 
R = Ra+Re- From Eq. ([H), one identifies R^ = - 2B^, - and R^ = -2A5 - 6AC. No- 
tice that by our approach we compute F^ —Ra as a single term, so that it is not possible to separate 
the two contributions. 

In other formulations , ^^'^^ the derivative of the energy with respect to the occupation number 
of a broad band is computed by shifting the band with respect to the Fermi level. This adds a non- 
interacting contribution to the curvature of the energy, since the level whose occupation is varied 
is itself a function of the band occupancy. Such non-interacting term has to be subtracted when 
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computing U by these approaches. ^ Our formulation is conceptually more similar to scaling the 
occupation of all valence atomic levels in a uniform way, and the non-interacting term is vanishing. 



C. Computational details 

The results presented in this paper have been obtained by DFT calculations with the Perdew- 
Burke-Emzerho^ generalized gradient approximation for the exchange and correlation func- 
tional. We used an all-electron linearized augmented-plane-wave code to perform the simulations 
with constrained core occupations. Periodically repeated supercells at the experimental lattice 
constants were adopted to describe the solids. One atom was ionized in a unit cell containing 
four and eight atoms for Cu and Zn, respectively. In both cases the ionized atom has no ionized 
nearest neighbours. Cell neutrality is preserved by increasing the number of the valence electrons, 
simulating the screening of the core hole by the solid. The spin-orbit splitting in core states as well 
as in the final state with two holes was taken into account by adopting DFT energy shifts for free 
atoms,^^ and is here assumed independent on the fractional charge q (we verified that the latter 
approximation affects our final results by no more than 0.2 eV). As for the coefficients A, B, and 
C in Eq. Q, we found values numerically more stable, with respect to convergence parameters, 
by performing a second order expansion of the eigenvalues rather than a third order expansion of 
the total energy. Therefore, we made use of Janak's theorem and Eq. ([7]), with eigenvalues relative 
to the Fermi level in the solid (hence, resulting XPS and Auger energies are given with respect to 
the same reference). Fulfillment of Janak's theorem and coincidence of results of Eq. ^ and (|7]) 
were numerically verified to high accuracy in a few selected cases. The values of q ranged from 
to 0.05 at intervals of 0.01. Comparison with denser and more extended meshes for the free 
atom case showed that results are not dependent on the chosen mesh. Matrix elements and Slater 



integrals and are taken from Ref. 



I2L and core hole lifetimes from Ref. |22|. 



III. RESULTS 



In this Section we report our results for the Auger lineshape of Cu and Zn. The 

core and the valence indices, c and v in the previous Section, are specialized to the 2p and 3d level 
of such elements, respectively. 

As an example of our procedure to extract the parameters A, B, and C [see Eq. Q], we report 



9 




0.02 0.04 0.02 0.04 

1 1 

FIG. 1: Example of fitting the Kohn-Sham eigenvalue to extract the coefficients A, B, and C. Results 
are shown for the 2p level of metal Cu. Panel (a) plots the KS eigenvalue (relative to the Fermi energy) 
with reversed sign (circles) and the fitted parabola from Eq. ^ (line), as function of the number of electrons 
removed from the 2p level, q. Panel (b) reports the same quantities after subtracting the linear term A + 25^. 

the case for the 2p level of Cu metal in Fig. [T] (the following considerations are also valid in the 
other cases). We remove the fractional number of electrons q from the 2p level of a Cu atom, and 
plot its (minus) KS 2p eigenvalue in Fig.[T^. Such a curve is fitted by the expression in Eq. (|7]). It 
is apparent from Fig. [1^ that a linear fit already reproduces the KS eigenvalue in this range of q to 
high accuracy. However, since results are to be extracted up to ^ = 1 or 2, the quadratic term in 
the expansion is also of interest. This is shown in Fig. [l}?, where the linear contribution (A + 25*5') 
has been subtracted. The parabola accurately fits the numerical results, with residuals of the order 
of 10 - 50 //eV. 

Table |n] collects our results for the coefficients A, B, and C, needed for the determination of the 
L23M45M45 lineshape of Cu and Zn. The values of B and C for the L2 and L3 cases are identical, 
following the assumption that spin-orbit splitting is independent of the fractional charge. The 
coefficients 5^45 and Cm^^ in the solid have been obtained by comparing results for core levels 
in the bulk and in the free neutral atom according to Eqs. (|10m II) . with the values of AS and 
AC averaged across the core levels as reported in Table HI Their negative sign indicates that the 
interaction amongst the two holes is more effectively screened in the solid. 

Let us now consider the dependence of our results on the particular choice of the reference 
atomic configuration. For comparsion, the Cu+ and Zn+ ions (with one electron removed from 
the As shell) have been used as a starting point for the evaluation of the atomic coefficients instead 
of the neutral one. We find similar modifications, canceling each other in Eq. (fTOl) . for core and 
valence B^ atomic coefficients (larger by about 1.5 eV in Cu and 1.3 eV in Zn). The same is 
found for the Cf coefficients (lower by 0.2 eV in Cu and 0.1 eV in Zn). As a consequence, the 
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Level 


A" B" C" 


A B C 


E^^^ Exp. 




'-'2 


930.17 27.74 1.25 


928.40 22.84 0.43 


951.67 952.0 


Cu 


L3 


909.81 27.74 1.25 


908.04 22.84 0.43 


931.31 932.2 




M45 


5.04 5.72 0.91 


2.86 0.87 0.10 


3.84 3.1 




L2 


1019.40 30.44 1.08 


1016.98 26.17 0.82 


1043.97 1044.0 


Zn 


L3 


995.69 30.44 1.08 


993.27 26.17 0.82 


1020.26 1020.9 




M45 


10.14 7.06 0.77 


7.53 2.82 0.42 


10.78 9.9 



TABLE II: Coefficients for the expansion of the total energy as a function of the number of electrons, 
E{N — q), for atomic {A", B", C) and bulk (A, B, C) Cu and Zn. As for the M45 values: by AM45 we indicate 
(minus) the weighted average of the 3d KS band; Bm^^ and Cm^^ obtained according to Eqs. dTOlITT]) . 
Theoretical XPS energies are given by Eq.®; experimental data are taken from Ref.ll2|. Values in eV. 

values for 5^45 differ by less than 0.2 eV, and those for CM45 are identical within 0.01 eV, with 
data reported in Table UIl Hence, as anticipated in the previous section, the choice of the reference 
atomic configuration does not affect significantly the evaluated XPS and Auger energies. 

Recall now that A + 5 + C is our estimate for the XPS excitation energies [see Eq. ([8])], which 
are reported in Table UIl and compared with experimental values.— Notice that bare KS excitations 
energies can be 30 eV smaller than the experimental value, but the addition of B and, to a smaller 
extent, of C, properly accounts for the missing relaxation energy, the left discrepancy being smaller 
that 1 eV. 

We report next our results for the 3d component of the IhLDOS, J((o) , for Cu and Zn in Fig. [2l 
We remind that such quantity is obtained by converting the KS density of states into the hole 
picture, and by translating the result by By + Q, to account for relaxation effects [see Eq. Q]. The 
total d IhLDOS, (i((o), is shown as a shaded area together with the components on the different 
irreducible representations over which the d matrix is diagonal. For both metals, the various 
components differ among themselves in the detailed energy dependence, but their extrema are 
very similar. 

As a final ingredient. Table Hill lists the values of U for the five LS components of the multiplet, 
computed by Eq. ([T4l) . Notice that the inclusion of a q dependence in U (via a cubic term in the 
expansion of the total energy with respect to a fractional charge) proves to be quite important. 
Indeed, in evaluating U the C coefficient is counted six times, hence bringing a larger contribution 
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10.5 11 

Energy (eV) 



FIG. 2: One-hole local density of states, (i((o), for Cu (top) and Zn (bottom), relative to the Fermi energy 
and normalized to unity. The shaded area is the total d DOS, d. 









Cu 


2.38 


7.76 3.34 2.67 2.25 0.33 


Zn 


8.16 


14.31 9.26 8.49 8.01 5.82 



TABLE III: Values of U resulting from the application of Eq. ([141 ). in eV. Slater's integrals from Ref. 
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than in the XPS energies previously discussed. Such inclusion gives an estimate ofU = U{q = I) 
which is 0.60 eV and 2.52 eV larger for Cu and Zn, respectively, than the corresponding values 
obtained as U{q = 0). 

We then compute the L23M45M45 Auger spectrum following Eq. dH). The outcome has been 
convoluted with a core hole lifetime of 0.49 and 0.27 eV (0.42 and 0.33 eV)2^ for the L2 and L3 
lines of Cu (Zn), respectively, and results in a multiplet of generally narrow atomic -like peaks, 
shown in Fig.|3l 

To analyze these results, let us focus on the principal peak (^G) in the spectrum, which can 
be associated with the absolute position of the multiplet. (The internal structure of the multiplet 
in our description only depends on the values of and F'^ which, as previously specified, were 
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FIG. 3: Simulated L23M45M45 spectrum for Cu (top) and Zn (bottom) metals. The vertical lines mark the 
position of the principal (^G) peaks from experiments.— 

taken from the literature.) The experimental energy of the (most intense) transition is marked 
by a vertical line in Fig. |3l The agreement of our results is rather good considering the absence of 
adjustable parameters in the theory: focusing on the L^VV line, the peak position (918.0 eV 
from experiments) is overestimated by 1.6 eV, while the one for Zn (991.5 eV) is underestimated 
by 1.9 eV. 

IV. DISCUSSION 

It is interesting to compare these results with those obtained by an expression commonly 
adopted for Auger energies, i.e., Qd^s ~ £c — 2£v — Uis- This is an excellent approximation when U 
is larger than W: for example in Zn, where W ^ 1.5 eV and Uiq = 9.26 eV, its application to the 
computed parameters yields a value which is only 0. 10 eV larger than the peak position derived 
from Eq. (|4]). However, when U is of order of W, significant deviations can be observed: e.g., for 
Cu (W ~ 3.5 eV and Uiq = 3.34 eV) the ^G position is overestimated by 0.66 eV. For smaller 
values of U, the quasi-atomic peak is lost for a broad band-like structure. This is the case for the 

component of Cu (the rightmost shoulder in the spectrum) for which we obtain Uip = 0.33 eV. 
However, this is an artifact of our underestimate of i7sph in Cu: experimentally, the peak is 
resolved as well. 

Besides these observations, the expression cO/,5 ^ £c — IZy — Uls is accurate enough for the 
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£c 


-2s,, 








Theory 


931.31 


-7.67 


-3.34 


920.29 


Cu 


Exp. 


932.2 


-6.2 


-8.0 


918.0 




Diff. 


-0.9 


-1.5 


4.7 


2.3 




Theory 


1020.26 


-21.55 


-9.26 


989.45 


Zn 


Exp. 


1020.9 


-19.8 


-9.5 


991.5 




Diff. 


-0.6 


-1.8 


0.2 


-2.1 



TABLE IV: Decomposition of the LT,Mi^^Mi,^ Auger kinetic energy into its contributions, according to 
the simple approximation co = — 28,, — U . Theoretical values of s^, £v. and U from tables [II] and |llll 
experimental data from Ref. 12. Values in eV. 



peak to discuss the discrepancy of our results with respect to the experimental ones. Let us focus 
on the L^yV part of the spectrum, and decompose the Auger kinetic energy co into its contributions 
(see Table HVl). Despite the fact that the overall agreement is similar in magnitude for Cu and Zn, 
it is important to remark that this finding has different origins. In both metals, we underestimate 
slightly the core photoemission energy and overestimate the valence photoemission energy by a 
similar amount. Both effects contribute underestimating the kinetic energy. In Zn, where our 
value of U is excellent, the error in (0 stems from the errors in the photoemission energies. In 
Cu, instead, U is seriously underestimated. This overcompensates the error in the photoemission 
energies, resulting in a fortuitous overall similar accuracy. 

We now examine the relative weight of two different ingredients of our method. First, the role 
of the spin-orbit interaction in the final state, which will be analyzed by comparing with results 
where such term is neglected; second, the resolution of the 2hLD0S in its angular components. 
To this respect, we notice that the matrix expression for the spectrum given in Eq. dH) can be 
significantly simplified under the assumption that the 2hLD0S is spherically symmetric and the 
spin-orbit contribution can be neglected. In this case, we can just take the spherically averaged, 
i.e., the total d IhLDOS, and compute its self-convolution, D^'^^ = d*d. An averaged Green's 
function, G^^\ is then defined as the Hilbert transform of D^^\ By replacing G^^^^j^^ l'S'j'm' ™ 
Eq. (lU with the diagonal matrix ^^sjMjMS' j'm' 

, we obtain the simple scalar equation 



5((0) = — y{2J+\)\ALs\hm 



LSJ 



(0 



l-t/L5G(o)(£,-a)) 



(15) 
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scalar Eq. (15) 
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Kinetic energy (eV) 

FIG. 4: Simulated spectrum for the L3M45M45 line of Zn. Solid line: full treatment of Eq. (01), as presented 
in this paper. Dashed line: neglecting the spin-orbit interaction in the two-hole final state. Dotted line: 
adopting the spherically averaged 2hLD0S and the scalar formulation, Eq. ( fT5] ). The origin of the vertical 
axis is shifted for improved clarity. 

where the dependence on LS quantum numbers is only via the matrix elements and the interaction 
matrix U, and each LS component of the spectrum is decoupled from the others. 

The Auger spectra simulated neglecting the spin-orbit interaction and calculated with the sim- 
plified expression of Eq. (fTSi) are plotted in Fig. |4]as dashed and dotted line, respectively, to be 
compared with the result of the full calculation [Eq. dH)], solid line. For simplicity, we limit the 
discussion to the L^VV line of Zn. The resemblance of the three results is remarkable. Indeed, in 
Cu and Zn the spin-orbit splitting for 3d levels is relatively small, 0.27 and 0.36 eV, respectively.^- 
Furthermore, despite the differences which characterize the angular components of the IhLDOS 
(see Fig. |2l), the convoluted 2hLD0S are only mildly different from D^^\ as reported in Fig. [51 
Now, in systems with a large U /W ratio, fine details of the 2hLD0S are not relevant for the po- 
sition of quasi-atomic peaks, which only depend on the weighted averages E^g as in Eq. ([5]). In 
our case, the values of E^g lie within 0. 1 eV from those corresponding to the averaged 2hLD0S. 
Hence the practically equivalent results obtained by Eq. (01) and Eq. (fT5l) . 

This analysis shows that, for a wide class of systems with strong hole-hole interaction, weak 
spin-orbit interaction, and a spherical symmetry to some extent, the simple formulation presented 
in Eq. (fTSi) is practically as accurate as the expression in Eq. (H)). One should instead adopt the 
full treatment for, e.g., heavier elements, or systems with low dimensionality. This remark is 
independent of the methodology to determine the parameters entering the model, either fully ab 
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FIG. 5: Non-interacting two-holes density of states, D^^\ for Cu (left) and Zn (right). The solid curve 
represents D^^^ = d*d, the self-convolution of the total d IhLDOS; the shaded area indicates the largest 
deviations from this result found amongst the angular resolved 2hLDOSs. 

initio as in the present approach, or by phenomenological arguments. 

Our method provides an agreement with experimental photoemission energies of the order of 
1 eV and of 2 eV for Auger energies. We consider this to be a rather good result, as a starting 
point, considering the absence of adjustable parameters in the model which is the new feature of 
our approach for CVV transitions in correlated systems. We trust that our simple method could 
already yield qualitative information on the variations of the spectrum to be expected following to 
modifications in the sample, e.g., when the emitting atom is located in different enviroimients. Of 
course, a much better agreement would be obtained by inserting phenomenological parameters, 
but at the cost of loosing predictive power. 

The internal structure of the lineshape, i.e., the multiplet splitting, is given very precisely. 
However, the first-principles treatment of the latter is not a new aspect of our approach, which 
is indeed based in this respect on atomic results available in the literature since decades. Let us 
instead focus again on the estimated position of the multiplet, which crucially depends on the 
parameters evaluated by our ab initio method. Even though the agreement with the experiment 
is about as good as for Zn as for Cu, actually the results for Zn are much better. In Zn, the 3d 
band is deep and narrow, and electronic states bear mostly atomic character. Their hybridization 
with the states closer to the Fermi level, which mainly contributes to screening, is small, somehow 
in an analogous way as for core states. Consequently, the approximations to neglect the energy 
dependence of By and Cy, and to transfer the values of AB and AC from the core states to the 
valence ones, produce very good results. In Cu, instead, the 3d band is higher and broader, and 
hybridizes significantly with the s-like wavefunctions. Our approximations turn out to be less 
adequate: the resulting U is about half the one derived from experiments. 
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Part of this discrepancy might also have a deeper physical origin, since the Hamiltonian 
adopted, Eq. ([T]), does not allow for the interaction amongst two holes located at different atomic 
sites. The CS model has been extended to consider the role of interatomic ("off-site") correlation 
effects which mainly produce an energy shift of the Auger line to lower energies.^ Studies based 
on phenomenological parameters suggest that such an energy shift could be of about 2.5 eV in 
Cu^ (smaller values are expected in Zn where holes are more localized and screening is more 
effective). For sake of simplicity, the off-site term has not been considered here and is left for 
future investigations. The parameters entering this term could be determined by ab initio meth- 
ods in analogy to the procedure shown here for evaluating U. It is however important to notice 
that adding the off-site term would not fix all the discrepancies observed in Cu, where also the 
lineshape, in addition to the peak position, is not satisfactory owing to small values for the on-site 
interaction U (e.g., the non-resolved peak). 

Enhancing the accuracy of the values of U seems therefore the most important improvement 
for the method presented here, especially for systems with broad valence bands. As a possibility, it 
would be interesting to use approaches which are capable to compute the total energy in presence 
of holes in the valence state. The methodology presented in Ref. 19, in which the valence occu- 
pation is changed by means of Lagrange multipliers associated with the KS eigenvalues, could 
be particularly effective. One should pay attention as some arbitrariness is anyway introduced. 
Namely, the value of U does depend on the chosen form of the valence wavefunctions. Such an 
arbitrariness is compensated in LDA-i-U calculations performed self-consistently.^^ Furthermore, 
to apply this method to systems with closed band lying well below the Fermi energy, large shifts of 
the KS eigenvalues would be needed to alter the occupation of the band to an appreciable amount. 

Another possible improvement concerns the photoemission energies. Calculations by the GW 
approach^^ of the IhLDOS could be used to account for relaxation energies, rather than adopting 
Eq. ([8]). Results available in the literature (e.g., for Cu^) are very promising in that sense. It is 
interesting to notice that the factor 5 + C plays the role of a self-energy expectation value, and that 
the use of a single value of 5 + C to shift rigidly the band is formally analogous to the "scissor 
operator" often introduced to avoid expensive self-energy calculations. The accuracy of such rigid 
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shifts for valence-band photoemission in Cu is discussed in Ref. 

Systems with larger band width or smaller hole-hole interaction would require to extend the 
approach to treat the dependence of the matrix elements on energy together with the interaction 
in the final state. Releasing the assumption that matrix elements equal the atomic ones, as in 
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the current treatment, or that particles are non-interacting in formulations accounting for such an 
energy dependence (like, e.g., the one in Ref. m, would allow switching continously between 
systems with band-like and atomic-like spectra. This possibility is currently under investigation. 

Finally, let us recall the basic assumption considered here that the valence shell is closed, which 
is crucial to the CS model in its original form. Efforts have been devoted towards releasing this 
assumption, resulting in a formulation by more complicated three-hole Green's functions,^^ for 
which no ab initio treatment is nowadays available to our knowledge. 



V. CONCLUSIONS 



We have presented an ab initio method for computing CVV Auger spectra for systems with 
filled valence bands, based on the Cini-Sawatzky model. Only standard DFT calculations are 
required, resulting in a very simple method which allows working out the spectrum with no ad- 
justable parameters. The accuracy on the absolute position of the Auger features is estimated to a 
few eV, as we have demonstrated by the analysis of Cu and Zn metals. We have shown that in these 
systems further simplifications like neglecting spin-orbit interaction for the two valence holes, or 
the non sphericity of the emitting atom, give results practically equivalent to the full treatment. 
Attention should be paid to the problematic parameter U. We obtained such a term with a good 
accuracy for the more localized, atomic-like valence bands in Zn, while it results underestimated 
in Cu. Its occupation number dependence, included via a cubic term in the expansion of the total 
energy, has been considered, and shown to play an important role. 

This step towards a first-principles description of CVV spectroscopy in closed-shell correlated 
systems enables identifying improvements which future investigations could focus on. In particu- 
lar, one would benefit from detailed calculations of the single-particle densities of states (e.g., by 
the GW method), from truly varying the valence-band occupation to obtain the U parameter, and 
from including off-site terms in the Hamiltonian. Prospectively, it would be desirable to take into 
account the energy dependence of the transition matrix elements. 
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